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Abstract 

We analyze a self-consistent theory of crack growth controlled by a cumulative damage variable d(t) dependent on 
stress history. As a function of the damage exponent m, which controls the rate of damage dd/dt oc a m as a function of 
local stress a, we find two regimes. For < m < 2, the model predicts a finite-time singularity. This retrieves previous 
results by Zobnin for m = 1 and by Bradley and Wu for < m < 2. To improve on this self-consistent theory which 
neglects the dependence of stress on damage, we apply the functional renormalization method of Yukalov and Gluzman 
and find that divergences are replaced by singularities with exponents in agreement with those found in acoustic emission 
experiments. For m > 2, the rupture dynamics is not defined without the introduction of a regularizing scheme. We 
investigate three regularization schemes involving respectively a saturation of damage, a minimum distance of approach 
to the crack tip and a fixed stress maximum. In the first and third schemes, the finite-time singularity is replaced by a crack 
dynamics defined for all times but which is controlled by either the existence of a microscopic scale at which the stress is 
regularized or by the maximum sustainable stress. In the second scheme, a finite-time singularity is again found. In the 
first two schemes within this regime m > 2, the theory has no continuous limit. 

1 Introduction 



c3 



The fracture of materials is a catastrophic phenomenon of considerable technological and scientific importance. Despite 
J> the large amount of experimental data and the considerable effort that has been undertaken by material scientists [p, many 
^ questions about fracture remain standing. There is no comprehensive understanding of rupture phenomena but only a partial 
classification in restricted and relatively simple situations. This lack of fundamental understanding is reflected in the absence 
of reliable prediction methods for rupture based on a suitable monitoring of the stressed system. 

Some progresses have been obtained in recent years in the Physics community. Based on analogies with phase transi- 
tions, several groups [p|]-[^2|] have proposed that, in heterogeneous materials with disorder such as fiber composites, rocks, 
concrete under compression and materials with large distributed residual stresses, rupture is a genuine critical point, i.e., the 
culmination of a self-organization of diffuse damage and micro-cracking characterized by power law signatures. Experi- 
ments f7[ |9, 11, |2|], numerical simulations [|], |6], [8, 10] and theory [|| confirm this concept. 



As a signature of criticality, acoutic emissions radiated during loading exhibit an acceleration of their rate close to 
rupture [Q, [J [12|]. Specifically, under a constant stress rate, the cumulative acoustic energy E{t) released up to time t can be 
expressed as 

E{t) =Eq- B{t c - t) a , (1) 



with B > and < a < 1. Expression ([]]) corresponds to a rate dE/dt of acoustic energy release diverging at the 
critical rupture time t c . This behavior ([]]) has been at the basis of previous claims that rupture is a critical phenomenon. 
In addition, this power law (|T]) as well as extensions with log-periodic corrections have been suggested to be useful for 
prediction [|13|, |, |, [TTJ, [I2p. 

Our purpose here is to present, extend and analyze a simple self-consistent model of damage that predicts a behavior 
similar to ([I]). We explore its different regimes and then improve on its "mean-field" version which predicts an unrealistic 
finite-time singularity. In this goal, we propose to use the general functional renormalization approach developed by Yukalov 
and Gluzman to cure this anomaly. We show how this technique allows us to change an unrealistic singularity into the 
observed behavior ([[]) with a reasonable exponent z = 1/2, without introduction of any extra parameters in the theory. 

2 Cumulative damage model 

Initially introduced as a global "mean field" (uniform) description of the global deterioration of the system at the macroscopic 



scale []14|], the concept of "damage" has been extended at the mesoscopic scale to describe the heterogeneity and spatial 
variability of damage in different locations within the material [[If], |l5[ [l7|]. We use the formulation of Zobnin Jl8| ] and 
Rabotnov fll9|] to show how it leads naturally to a finite-time singularity. We first recall briefly the integral formulation of 



Rabotnov [ 19] (pages 166-170) and then transform it in differential form to exhibit the fundamentally nonlinear geometrical 
origin of the singularity. 

A material is subjected to a stress o"o at large scale and each point r within it carries a damage variable d(r, t). When 
d reaches the threshold d* as some location, this local domain is no more able to sustain stress and a microcrack appears, 
leading to a redistribution of the stress field around it according to the laws of elasticity. The local damage d(r, t) at point r 
at time t is supposed to evolve in time according to 

^ = k(r,f)P, (2) 
where a(r, t) is the local stress field at point r at time t and m is a damage exponent which can span values from to 



close to oo depending upon the material. In the discrete 2D models of Refs.[g, |15[ ], it was shown that rupture reduces to the 
percolation model in the limit m — > 0. In the other limit, m — > oo, rupture occurs through a one crack mechanism. 

Following Rabotnov [|TJ], we assume that a major crack dominates the rupture process. If only one crack is present 
within the system, the stress a(r, t) is easily calculated. Considering only the possibility of a linear straight crack of half- 



length a(t) advancing within the material at velocity da/dt (see [20] for generalizations to self-affine crack geometries), it is 
enough to calculate the stress field on the points ahead of the crack to fully characterize the rupture dynamics. For a planar 
elastic material subjected to a uniformly distributed antiplane stress at infinity with a crack lying on the y-axis between 
—a(t) and a(t), the stress field at point y on the y-axis beyond the crack tip is 

= ~3~ /P - \a(t)\' ' 0) 

The mean-field approximation made in this first version of the model consists in assuming that the stress field is not mod- 
ified by the non-vanishing and non-uniform damage field. This means that the elastic coefficients are taken constant and 
independent of the progressive damage, except of course when the damage reaches its rupture threshold d* . 

The law describing the growth of the crack, i.e., the dynamics a{t), is obtained from the self-consistent condition that 
the time it takes from a point at y, at the distance y — a{r) from the crack tip at time r, for its damage to reach the rupture 
threshold d* is exactly equal to the time taken for the crack to grow from size a(r) to the size a(t) = y so that its tip reaches 
the point y exactly when it ruptures. This is illustrated in figure 1 . Mathematically, this self-consistent condition is that the 
integral of (HJ) from time at which the pre-existing damage was till time t at which the crack tip passes through y is such 
that d reaches exactly the threshold d* at the time t. Two conditions must thus be verified simultaneously: 

1. a(t) = y (the crack tip reaches point y) and 

2. d(y, t) = d* (the damage at y reaches the rupture threshold). 



3 The linear damage law: m = 1 



We first consider the linear damage law m = 1 corresponding to the initial formulation of Zobnin [|18|]. This case has also 
been investigated and solved in [21] in the context of crack growth due to electromigration rather than mechanical stress 
(the current plays the role of the stress and, in the antiplane case studied here, the two problems are formally identical). 
This model is particularly interestingsince it allows both for an exact solution and an exact renormalization in the functional 
renormalization scheme [30]. It also provides a benchmark for approximate solutions in the general case < m < 2 as we 
discuss below. 

We now proceed to give the equation for the crack dynamics and its solution. By integration of (||), the two self-consistent 
conditions expressed for the case m = 1 lead to 

- f> - = d% (4) 

where the loading stress <7q can depend on time. The solution of this integral equation provides the time evolution a(t) of 
the macro-crack. To get it explicitely, we set 

z = [a{t)} 2 and ( = [<*(t)] 2 • (5) 
Changing the variable of integration from r to ( gives 

, r (dr/dQ 3d* 

«C °o , 7 = T.—j= ■ (o) 

This equation (Q) is an Abel equation with index —1/2, involving a fractional integral operator Defining the Abel 
operator I* acting on the function f(t) as 

Iaif} = f ^TX^ f( s )ds , (7) 
Jo r(l + a) 

the product of two such Abel operators is 
This shows that 

We thus see that /* J* i_ a {f} = Iq, which is nothing but the integral operator. The inverse of the Abel operator I* is thus 
(J*i_ a ). Applying this result to (H), we find 

dr _ 3 d* d f* d( 



dz 2 7r<7 dz J Zo y/((z - C) 

Calculating the integral in the r.h.s. of (JlOt), performing the derivative and inverting to get dz/dr, we get 

dz 2tt(Tq Z\Jz — zq 
dr 3d* ■ s /~Zq 



(ID 



Replacing z by [a(i)] 2 leads to the differential equation for the crack half-length a(t) 




which is exactly equivalent to the self-consistent integral equation @. It is remarkable that the local growth equation ([H 
embodies exactly the same physics as the long-term memory integral (^). 

For simplicity, let us take the loading stress <to constant. This situation is generic of experiments measuring the lifetime 
of structures under a constant load. At sufficiently long times for which a(t) » oq, expression (|l2|) reduces to 

da vro-o 2 ni\ 
— ~ — a . (13) 

at 3ao«* 



Equation (13) is characteristic of a solution going to infinity in finite time. Indeed, we can write (13) as da/dt oc ra, with 
a growth rate r oc a. The generic consequence of a power law acceleration in the growth rate r oc a s with 5 > is the 
appearance of a singularity in finite time: 

a(t) oc (t c - t)' 13 , with./3=- and t close to t c . (14) 



Equation ( |13| ) is said to have a "spontaneous" or "movable" singularity at the critical time t c Q23Q, the critical time t c being 
determined by the constant of integration, i.e., the initial condition a{t = 0) = ao- Note the intriguing fact that the (t c — i)" 1 
singularity appears as the solution of a linear mechanical problem. The source of the quadratic nonlinearity is the non-local 
geometrical condition that the delayed action of the stress field on the cumulative damage should coincide exactly with the 
passage of the crack tip. The nonlinear finite-time singularity has thus fundamentally a non-local geometrical origin, or 
alternatively can be seen to result from a long-term memory effect. 
The exact solution of ( |I2| ) is easily obtained by integration: 

This retrieves the solution obtained by Zobnin [18] and Rabotnov Jl9|]. We verify directly that the singularity occurs when 
the cosine goes to zero, i.e., when the argument reaches tt/2, i.e., for t c = 3d* /2<jq. Since the cosine vanishes linearly with 
time, this recovers the asymptotics ( |I4| ) with the exponent j3 = — 1, as predicted by the asymptotic equation (|T3|). 

4 The nonlinear damage law with finite-time singularity: < m < 2 
4.1 Derivation of the differential equation for the crack dynamics 



The case where < m < 2 can be similarly treated and our results here extend those of Zobnin [hq] and Rabotnov [19]. Our 



results retrieve those found in [24], obtained in the context of crack growth due to electromigration. For completeness and 
coherence in notation, we briefly present the method and the results which are more focused on the finite-time singularity. 
For simplicity, we impose <ro constant. Integrating (||) and applying the self-consistent conditions 1-2 leads to 

* (2a \ m [a(t)] m 



dr — ^ ^ sr = d* . (16) 

o v 3 ; ( [a (t)]2 _ [o(t)]2)t 

We set again the change of variables (Q) and changing the variable of integration from r to £ gives 

rfc ^U(^Y"4. (17 , 



This equation @ is again an Abel equation with index — m/2 if < m < 2. 

In order to transform it into differential form, we could use the formalism of Abel operators. We choose a more trans- 
parent and direct approach which is closely related. First, we multiply both sides of (|T^) by l/(y — z) l ~~ and integrate 
over z from zq to y: 



dz 



20 



d( 



(dr/dC) 



(z-C) 2 (y-z) 1 a 
Changing the order of integration in the l.h.s. leads to 



3 
2^ 



d* 



dz 



z 2 (y — z) 2 



(18) 



z d( 



dz 



C {Z-C)2 (y-z) 



I rn 

1 2 



3 
2^ 



d* 



dz 

m ' \ 1 m ~ 

z z 2 (y - z) L 2 



(19) 



where we have used the equality of the triangle fj£ dz J5 = J 2 2 Q c?C dz. 

The integral in the square bracket in the l.h.s. of ( |l9| ) can be expressed through the Euler beta-function B(a, b): 



dz 



We thus obtain 



C (z-C) 2 (y-zY 



7T 



B 1 



r(i 



t) r(f) 



7T 



r(i 



sin (mf ) 

After differentiation with respect to z, we get 



2cr r 



m _i_ mi 
2 ' 2 / 



sm m4 



rf2 



Gao) 



m sin (m|) d 
7T <iz 



m . -i m 

Zo Z2 (y - z) 1 2 

C -( Z _ C )1-- 



(20) 



(21) 



(22) 



4.2 Asymptotic solution close to the finite-time singularity 

The solution of (|22|) can be obtained for large crack sizes a(t), i.e., large z. In this goal, we replace the term (z — C) 1 T 
in the integral in the r.h.s. of q22|) by z a , neglecting £ compared to z. Intuitively, this is justified over the whole domain 
of integration because the contribution from the domain where ( is not negligible compared to z is finite, since the power 
1 — Y is less than one, corresponding to an integrable singularity. 

With this approximation, the integral can be performed, the derivative taken and after inverting, we get 



\2a 

Using a = ifz as defined in (||), we obtain 



d* 



sm [m 



f) dz 



7T 



where 



and 



a(t) 



/? = 



d,T 



a 



Z Q 2 Z Z 2 



(1-4)' 



m , 2 sin (m|) 



vr(2 - m) 



(23) 



(24) 



(25) 



(26) 



2 — m 

Note that the exact asymptotics ( [14] ) of the case m = 1 previously solved exactly is recovered, with the correct exponent 
z(m = 1) = 1 and a rather good approximation of the critical t c : while the exact value is t c = 3d*/3o"o, expression (25) 



predicts j-t c , i.e., 36% lower. The critical time t c as a function of m is smooth with no accident or divergence over the whole 
interval. In particular, the estimated critical time for the limit m — > 2~ is equal to 9d* /4<7q. 

In contrast, the exponent z increases from z(m — ► + ) = 1/2 to +00 as the damage exponent m varies from to 
2. The limit z(m — ► + ) = 1/2 can be rationalized as follows. This limit m — > + corresponds to the situation where 
damage becomes independent of stress. As a consequence, reintroducing some heterogeneity for instance on the pre-existing 
damage, rupture is then equivalent to percolation, as the parts of the system that break as a function of time are determined by 
the damage accumulating at the same rate for all point but with different random initial values. In mean field percolation [25] 
obtained through the consideration of one-dimensional percolating paths consistent with the present one-crack geometry, the 
elastic energy under constant load diverges a (t c — i) -1 where 1 is the mean field value of the exponent t for conductivity 
(which is the same as elasticity in the scalar mode III version of mechanical deformations used here). Since the elastic 
energy is proportional to the square of the crack length, we get the prediction a(t) ~ (t c — t)^ 1 / 2 . This reasoning holds if 
the exponent is a smooth function of disorder and geometry (the present studied here is a the zero-disorder limit). 

The divergence of z at m = 2 signals a change of regime that we study in the next section. 



5 The nonlinear damage law with m > 2 



For m > 2, the integrals in the equations ( |i6| ) and ( fr7j ) diverge at a(r) = a(t), since the negative power with exponent m/2 
is no more integrable. Technically, the main difference between the cases m < 2 and m > 2 is that we need to regularize 
the infinity in the expression (^0|) by introducing some sort of dimensionless cut-off. The important physical message is that 
the regime where m > 2 is controlled by a novel physical parameter, which we identify as a length scale associated with 
the damage law. In other words, the physics of the rupture is inherently controlled by the choice of the cut-off, i.e., by the 
existence of a microscopic length scale. We could summarize the situation by saying that there is no continuous limit to the 



theory for m > 2. This is similar to previous observations obtained in a dynamical theory of rupture front propagation [26]. 
We now present two ways for regularizing the divergence and thus for obtaining a meaningful theory of rupture. 



5.1 Regularization by damage saturation at a microscopic scale 

Before describing the physical content of the regularization we propose, we need to express the problem in a more manage- 



able mathematical form. Since the culprit for the divergence is the integral (g0() and the divergence occurs for ( — » z , we 
introduce the variable 

(27) 



Z 



y-C 



and rewrite ( 2C ) as 



dz 



c (z-c) 2 {y-z) 1 2 



dZ Z~~ (1 - Z) 



(28) 



which makes apparent that the divergence is due to Z 2 at the lower bound 0. It is thus natural to regularize by introducing 
a dimensionless cut-off e > and replace ( ^8| ) by 



dZ Z~~ (1 - Z)- 



Tft 1 



The function b(m, e) is such that 



lim t 



+0+ 



b(m,e) = B(l 



m m s 



b(m, e) 



for < m < 2 , 



where the beta function B(l — ^) has been defined in (20) 



In constrast, we have 



b(m, e) 



m — 2 ' 

e 2 



for m > 2 , 



(29) 



(30) 



(3D 



and 



b(m, e) ~ In - , for m = 2 , 



(32) 



showing that the divergence of the integral ( |28| ) is now encapsulated in the dependence of the factor b(m, e) on e. This 
regularization scheme thus relies on the existence of the definite integral (^), by analogy to the case m < 2. 
Using the regularization (|^), we obtain 



dr 



d* d 



dz \2a J b(m,e) dz J zo (~ (y - C) 1 



" 2 



which extends (10) to the regime m > 2. Its formal solution obtained in implicit form is 

d* 



t 



Jz„ dy J z . 



(JL 

V2cr y b(m,e) J Zo u dy J Zo (~ ( y - Q 1 ~ 



(33) 



(34) 



This regularization scheme allows to obtain exact solutions for integer m's. We examine the solutions for m = 2, 3, 4 
and 5 and then the general case. For m=2, we have the expression for all times given by 



where 



z 2 (t) = zoe 1 ' 1 ^ , 
9d* 1 



t. 



e,2 



4a 2 Ini 



For m = 3, z^(t) is the solution to 



\J(z/zq)-1 - taa-\y/(z/zo)-l) 



t 



where 



For large times, we get 



27d* 



^,3 



8<7 3 



£2 



For m = 4, z^ (t) is the solution to 



ln(z/z ) + (z/z ) - 1 



t 



where 



For large times, we get 



t 



e,4 



81ri* 
16cj 4 



04 (t) « Z — • 
te,4 

For m = 5, is the solution to 

((s/*o) - 1) 3/2 - 3 ((z/* ) - 1) 1/2 + 3tan- a [((z/z ) - 1) 1/2 

where 
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*e,5 



243d* 

32ct 5 



c 3/2 



(35) 
(36) 

(37) 
(38) 

(39) 

(40) 

(41) 
(42) 

(43) 
(44) 



For large times, we get 



More generally, at large times 




(45) 



2 

n-2 



z m (t) f»Zo[- — , (46) 

( t \ 

a m (t) w ao - — , (47) 



where 



m-2 1 

t e m oc e 2 ~ _ . (48) 

b[m, e) 

From these solutions, it is apparent that the dynamics z(i) = [a(t)] 2 is controlled by the characteristic time t e>m defined 



in (|48j). Note that the inverse dependence of t e ^ m on b(m, e) is obvious from the expression (J34Q. As the cut-off e — > 0, 
te,m —> and the global rupture occurs in vanishing time. The physical explanation of this phenomenon is as follows. For 
m > 2, the driving force a m of the damage law (Q) is so strong close to and at the crack tip, that it takes effectively zero 
time for a point to be brought to the damage threshold. To see this, let us truncate the integral in ( JT^ ) such that the upper 
bound is changed from ttot — r). The divergence of the integral at the crack tip means that the contribution to the cumulative 
damage occurring in the time interval from t — 7] to t is larger (actually infinitely larger) than the contribution from time to 
time t — rj. This means that the progressive damage leading to the acceleration of a(t) for m < 2 is replaced by an infinite 
velocity as soon as we start from a finite crack and do not introduce the finite cut-off length. 



This clarifies the physical meaning of the cut-off e defined in (|29|). A non-zero e means that the integral over £ in (|28|) 
does not go all the way up to z. Translated in terms of physical distances, it means that the integral in ( in the l.h.s. of 
(|T7|) also does not go all the way up to z. Physically, this means that the damage on a given point ahead of the crack tip 
reaches the critical value d* before the crack tip reaches that point. The value d* is no more the rupture threshold but a 
saturation value. The crack tip dynamics is now determined by the condition that the damage at any given point y reaches 
this saturation value d* when the crack tip is at a fixed distance oc e from y. This condition embodies the existence of a 
microscopic length scale oc e such that the damage is no more defined as smaller scales. 

The main result of our analysis is that the characteristic time scale t e>m of the crack dynamics is controlled by the 
microscopic length scale. The theory has thus fundamentally no continuous limit. It is one of several interesting and 
important examples in physics where the macroscopic physics is completely controlled by the microscopic physics (the ultra- 
violet cut-off). This situation is found in many physical problems, for instance in correlation functions in two-dimensional 



systems [|27|], in non-linear diffusion [ |29| ] as well as in quantum electrodynamics [[28|], Note however the difference between 
the last two examples and the former ones: in our rupture problem as well as in the case of correlation functions in 2D 
systems, the ultraviolet cut-off appears naturally, as an atomic distance, while in the last two case, there is no meaningful 



natural cut-off, hence necessity to "cover-up" divergencies by the "renormalization" procedure [28] 



5.2 Regularization by stress saturation at a microscopic scale 

The previous regularization scheme invokes a saturation of the damage at a microscopic length oc e. Alternatively, the 
saturation can occur on the stress field, whose mathematical divergence is bound to be rounded off at atomic scales. This 
provides another regularization scheme. To implement it, we use the continuous expression ^ for all distances from the 
crack tip down to a regularization length £ such that, for distances from the crack tip from to £, the stress is constant equal 
to a(£) given by (Q) with y = a(t) + £. This regularization is standard in the theory of damage and of plasticity. The idea 
is that a sufficiently large damage exerts a feedback on the stress field which then departs from its damage-free continuous 
expression (|3j> . This extension to Rabotnov's treatment provides a natural way for constructing a self-consistent theory of 
damage: not only does rupture occur by the cumulative effect of damage, damage has also the effect of smoothing out the 



mathematical singularity at the crack tip. The cut off £ has the physical meaning of a so-called process zone or damage zone 
and its introduction is fully consistent with the dynamical damage law However as we have seen above, for a damage 
law with an exponent < m < 2, the resulting dynamics becomes insensitive to the existence of a microscopic length scale 
in the limit where it is small. In this sense, the regime < m < 2 is more universal and has a continuous limit. 
We propose two models that implement these ideas. 



5.2.1 Saturation of the stress at a fixed distance to the crack tip 



The regularization scheme used here is such that the stress is assumed to saturate at a fixed distance £ from the crack tip. 
Thus, for times r up to tc(t), the damage at a fixed point that will be reached by the crack tip at time t is growing under the 
influence of the stress field created by the crack. From time tg(t) up to time t, the damage is increasing linearly with time, 
since the stress is assumed constant and equal to the value it reaches at time ti(t). The saturation time t^(t) is determined by 
the equation 

a(t) - a(t e ) = £ . (49) 



In this version of the regularized theory, expression ( |16| ) is changed into 

[a(t)r 



dT 



([a(tW ~ [a(r)f 



+ 



[a(t)Y 



[(2a(t) -£)£Y 



it - t e ) = d* 



(50) 



where ti is the time at which the crack tip is at the distance £ from the position it will have at time t (see (4-9)). 



Note that we now have two equations for two unknown a{t) and a{ti). The second term in the l.h.s. of ( pQ ) expresses 
the linear increase in damage from ti to t under the saturated stress. In this model, £ is fixed and t% ajusts itself. The value 
of the saturated stress is not a constant but increases as the crack gets larger and larger, since it corresponds to the value at a 
fixed distance £ from the tip of a growing crack. 

With the change of variables (0) and changing the variable of integration from r to ^ gives 



i^-ef (dr/dQ 



d£ 
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[z-Q- [(2^-£)£]T 
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a(tt 



Note that t — te can also be written 



t(z) -r((^-£) 2 ) « 2£j~z 



dr 
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(51) 



(52) 



The integral in the l.h.s. of ( |5l| ) is analyzed similarly to the previous case (|17|). We multiply the integral by l/(y - 



and integrate over z from (^/zq + £) to y: 
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J ■ 



dz 
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d C [J(^C+e) 2 {z-Q)' (y-zY 



(53) 



where we have used the equality of the triangle jf 



The integral in the bracket in the r.h.s. is the same as in (E9b, which defines the function b(m, e) with 



dzJ^dC 



(VC + £) 2 -( „ 2£V( 



Note that e(Q is now a function of 

Using ( |52| ) and (p3l), expression (|1|) gives 
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(55) 



Since 

m ^ 



b(m,e(Q) ocl/e^o, (^j) * , (56) 



we see that the first integral of the l.h.s. of ( p5| ) is negligible compared to the second integral of the l.h.s. of (55) in the limit 
of large cracks, i.e. large z. Neglecting i compared to sfz in the denominator of the integrant of the second integral of (55) 
and equating this second integral to the r.h.s. gives the following equation 



For m>2, ^ + i>l and the solution of (B7b is 



a(t) oc — , (58) 

(tc-t)' 



1 m — 2 



where t c is determined from the initial size of the crack. The finite-time singularity results from the ever-increasing stress 



field at the fixed distance I from the crack tip. This solution (^) is qualitatively different from the solution (£4J) found in the 
regime m < 2 as (|58|) depends in a fundamental way upon the existence of the regularization scale i. 

5.2.2 Saturation by fixing an absolute maximum stress 

An alternative prescription for the regularization is that the stress saturates at a constant value <7 max . This is in constrast 
with the previous regularization scheme where the stress saturates at a value reached at a constant distance, this value thus 
increasing with the crack length. Expression ( |50| ) is then changed into 

(59) 



^ /2a \ m [a(t)p 
Jo \ a J {[a(t)] 2 - [a(r)] 2 ) 2 



with 

2a a(t) 

-T, f = CT max , (6U) 

which is the condition that the stress saturates. It gives 

[a(t e )} 2 = [a(t)} 2 (1 - A 2 ) , (61) 

where 

A=^. (62) 

■J Cm ax 

The equation (|^) governing the dynamics of the crack tip can thus be written 

— d( J^tf + [<W] A * Tz = d , (63) 



where we have used the expansion 



t(z) - r((l - A 2 )z) ^A 2 z^\ z + (A*) , (64) 



valid in the interesting regime <r max » ctq giving A « 1. 



The second term in the l.h.s. of (63) dominates the first integral for large z, as can be checked a posteriori. For large 

nto 

[<x ma x] m A 2 z^- = d* , (65) 



crack sizes, the expression ( p4| ) can thus be simplified into 

dr 

dz 

whose solution is 

a(t) = a e t/t0 , (66) 

where 

2d* _ 9d* 

~ A 2 \rr l m ~~ In 2 \rr l m " 2 ' ( "* 

[°maxj zo o l°maxj 

6 Beyond the mean field version by functional renormalization 

Let us restrict our discussion to the case m = 1 for which we have the complete analytical solution for the crack dynamics. 
The solution (|l5|) with its asymptotic behavior ( |14] ) is not physically reasonable, as the crack reaches an infinite length 
in a finite time. The {t c — t)^ 1 singularity has been found to appear as the consequence of a geometric nonlinearity on 
an otherwise linearized mechanical problem. In reality, nonlinearity, viscosity, feedback, spatial heterogeneity of material 
properties and of cracking should modify the singularity. In addition, the main simplification in the previous approach is to 
neglect the impact of damage on the elastic coefficients of the material, thus leading to a stress field created by the crack 
which is identical to the field that the same static crack would generate in an undamaged material. Our hypothesis is that 
such modification can be deduced by a smooth or regular deformation of the solution previously obtained. 

In this goal, we propose to apply the Yukalov-Gluzman functional renormalization method [|3(]] to the series expansion 



of the solution ( |15| ) to obtain the renormalized law that accounts for these effects in a generic sense. Let us first consider the 
asymptotic power law singularity (|l4|) 

a(t) = — (1 - x)' 1 , where x = t/t c . (68) 

The powers x n in the expansion 

a(t) = — (l + x + x 2 + x 3 + ...) (69) 

may be considered as hidden free parameters. Indeed, let us multiply the expansion by x s . We then have a trial expansion 
for the solution. For s = 0, we return to the regular expansion. Such multiplication can be applied repeatedly, for instance 
using the functional renormalization method [[30|]. The idea behind the introduction of the multiplicative (control) function 
such as the power s in x s is to deform smoothly the initial functional space of the expression a(t) taken as an approximation 
to be improved. The condition for the improvement is to obtain a faster and better controlled convergence in the space of 
the modified functions upon addition of successive terms x n in the expansion. By this procedure, the dominant poles are 
eliminated or weakened as a result of a sequential reduction of stress level at each step of the resummation procedure. This 
corresponds to utilizing the information from the initial series pertaining to the times preceding the critical time t c , where 
the level of damage is lower. Thus, the renormalization procedure is performing a mapping from the dynamics at early time 
far from the critical point to later times closer to the critical time. The stabilization stems from the fact that the information 
contained in the initial series related to times close to t c is minimized on the basis that it has an overly destabilizing effect in 
the description and should be weighted less than the information at earlier times. 

At each step of the functional renormalization corresponding to the addition of a new term, we select the renormalized 
function according to the principle of minimum "local" multiplier, i.e., maximum stability on each sub-step of the renormal- 
ization procedure. Since these multipliers are proportional to the derivative da/dt [[3C|], the principle of minimal multiplier 
implies a selection of the real-time trajectory of the crack with minimal rate of damage (minimal stresses). In other words, 
this procedure amounts to improve the theory by allowing the crack to organize and develop so as to choose the most favor- 
able path or dynamics. It can be shown [|(]] that, at each step, the choice of a formally infinite exponent s corresponds to 
the minimal multiplier at arbitrary time. 



The functional form of a super-exponential solution is selected by this procedure: starting from an expansion 1 + a±x + 
CL2X 2 + (I3X 3 + ... + a k x , the renormalized expression is as follows. With the use of the notation 



bo = a , b h = —^—, k = 1,2,..., (70) 

we obtain the bootstrap self-similar approximant up to order k 

F k (x) = b Q exp(b 1 xexp(b 2 xexp(...b k ^ixexp(b k x)))...) , (71) 



introduced by Yukalov and Gluzman Q30|]. 

Let us now apply this result to the case (69) where all the coefficients a n are equal to 1. The corresponding renormalized 
approximant replacing the initial input 1/(1 — x) of the expansion reads 

F(x) = exp (xexp (x... exp (x)) ...) . (72) 

This embedded exponential series converges to a well-defined function. To determine it, we note that F(x) obeys the 
recursion relation 

F k+1 (x) = exp(xF k (x)) . (73) 
The fixed point to which these series of approximants converge is thus solution of 

F = exp (xF) , (74) 

noting F = na/2ao. The limit -Foo(x) exists for — e < x < e. 

The fixed point F{x) can be shown [^] to be the solution of the equation 

dF F 2 

Searching for a solution in the form of a Taylor series 

00 

F(x) = J2vnX n , yo = l, (76) 

n=0 

we get 

Vn = : • (77) 

nl 



Since n! n n e~ n , y n e n for large n and the generic term y n x n in the series ( pq ) is proportional to (ex) n . This shows 
that the radius of convergence of the series ( |76| ) is 1/e. 

F(x) has a singularity when x approaches 1/e from below, whose shape is obtained by expansions of expression (f74"l): 



F(x) = s _ 1/e e(l-y/2 e^^l/e-x^j 



(78) 



Thus, the self-similar functional renormalization has transformed a pole (divergence of o) at t = t c into a square root 
singularity (finite a) at a smaller t = t c /e. In this renormalized theory, the crack accelerates up to the time t c /e as which 
time its velocity diverges, while the crack is still finite. This announces the global breakdown. It is interesting that the 
exponent 1/2 is close to the value found for acoustic emissions in experiments ||, |TT], 12]. 

We can offer the following physical intuition for this transformation from the solution (14) with /3 = 1 to to f3 = — 1/2. 
As the material becomes more and more damaged, the ulterior functional dependence of damage as a function of applied 
stress is modified. Actually, the series of functional renormalization amounts to effectively evolve or renormalize the damage 
law (0) into a succession of effective laws captured by the sequence of approximants, each approximant order corresponding 



to an increase in the overall damage of the material. Here, we have a mapping between a measure of evolution via the 
cumulative damage, i.e., a measure of passed time, and the order of the approximants and thus the distance to the fixed point 
in the functional space. 

Consider now the general case a(t) ~ (1 — x)~@ . Expanding in power series, we get 



(1 - x)-P = J2 anx 11 , (79) 

n=0 

where 

an = n!(/?-l)! • (80) 
The Yukalov-Gluzman renormalization scheme gives the superexponential ([7l]) with coefficients b n given by 

b n = (n + P-l)/n. (81) 

Since b n — > 1 for large n for any [5, the fixed point of the approximants is controlled by the same finite square-root singularity 
of the type ([78]). Thus, the functional renormalization maps all finite-time singularities with different exponent f3 on the same 
universal law a(t) = a(t c ) — C^Jt c — t, where C is a constant depending in particular on /?. 



7 Concluding remarks 

Two main regimes have been found for the growth of a crack in a medium obeying the damage law d(d)/dt = a m (equation 
(§)), where a is the local stress. For < m < 2, a pre-existing crack grows to infinity in finite time and the divergence 
occurs as a power law finite-time singularity. For m > 2, the solution exists for all times but the characteristic time scale 
of the crack growth is an increasing function of a microscopic length scale, which is essential for regularizing the otherwise 
ill-defined problem. This microscopic length scale embodies the physical mechanism(s) by which the mathematical stress 
singularity at the crack tip of a perfectly sharp crack is rounded-off . We have examined two main scenarios, a damage-limited 
rupture and a stress-limited rupture. 

The remarkable behavior of this simple model results from the form of the irreversible damage law, in particular from 
the fact that any non-vanishing stress increases the damage. Damage at any point is thus a kind of ever increasing counter 
of the history of the stress on that point. This feature prevents the existence of stationary solutions of cracks propagating at 
constant velocities. In contrast, we only obtain "run-aways." 

Stationary solutions can be obtained in simple generalizations of the damage law (||), for instance with a stress threshold 
below which no damage occurs or with a healing or work-hardening term allowing recovering of the material and decrease 



of the damage when the stress is low. Such situations have been investigated in discrete two-dimensional models [15] 
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Figure 1: Illustration of the law governing the growth of the crack: the dynamics of its length a(t) is obtained from the 
self-consistent condition that the time it takes from a point at y, at the distance y — a(r) from the crack tip at time r, for its 
damage to reach the rupture threshold d* is exactly equal to the time taken for the crack to grow from size a(r) to the size 
a(t) = y so that its tip reaches the point y exactly when it ruptures. 



